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SUMMARY 

Static postbuckling and nonlinear dynamic analyses of plates are usually 
accomplished by multimode analyses, although the methods are usually com- 
plicated and do not give straightforward understanding of the nonlinear 
behavior. Assuming single-mode transverse displacement, a simple formula is 
derived in this paper for the transverse load-displacement relationship of a 
plate under in-plane compression. The formula is used to derive a simple 
analytical expression for the static postbuckling displacement and nonlinear 
dynamic responses of postbuckled plates under sinusoidal or random excita- 
tion. Regions with softening and hardening spring behavior are identified. 
Also, the highly nonlinear motion of snap-through and its effects on the over- 
all dynamic response can be easily interpreted using the single-mode formula. 
The theoretical results are compared with experimental results obtained using 
a buckled aluminum panel, using discrete frequency and broadband point excita- 
tion. Some important effects of the snap-through motion on the dynamic 
response of the postbuckled plates are found. 

INTRODUCTION 

If a plate is curved, initially or subsequently in service due to post- 
buckling stresses, the static and dynamic behavior in the transverse direction 
can be highly nonlinear, which may include hardening-spring, softening-spring 
or even an instability condition with snap-through motion. Theoretical and 



experimental results of large amplitude vibration of postbuckled plates under 
sinosoidal excitation were obtained by Yamaki and Chiba, 1 however, snap- 
through motion was not studied. The characteristics of snap-through motion in 
a postbuckled beam under sinusoidal excitation was studied by Tseng and 
Dugundgi. 2 A theoretical study of the random response of an initially curved 

O 

beam including snap-through motion was done by Seide. However, a thorough 
and straightforward understanding of the nonlinear behavior (particularly 
snap-through motion) of general curved plates is difficult to gather from the 
previous research results. The present study was conducted to fill this gap 
using a single-mode analysis method and experimental investigation with 
sinusoidal and random excitation forces on a postbuckled plate. 

GENERAL FORMULAE FOR NONLINEAR BEHAVIOR OF PLATES 
Equation for Equilibrium in the Transverse Direction 

For a plate under uniaxial compression with uniform edge displacement, 
the relationship between modal displacement and modal force for the buckling 
mode is given by (derived in Appendix I) : 

3 

for static equilibrium: q - Rq * p (la) 

for dynamic motion: q/u + 2tq/8 + (q - Rq) ■ p (lb) 

where 

q, nondimen8ional displacement parameter, * Q/Qp 
Q, modal displacement 
Qp, value of Q at R * 1 
R - A - 1 
A * u/u c 

u * in-plane edge shortening displacement 
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u c ■ value of u at which buckling starts 

ft * linear natural circular frequency of the flat configuration 
C “ modal damping coefficient 


p - nondimens lonal force parameter, - 

P - externally applied modal force 

K “ linear modal stiffness of the flat plate 

Qp , u c , ft, t , K depend on the assumed shape function of the mode 
and other plate parameters. The nondimensional parameters, q, R, P, can be 

evaluated after Q p , u c , K are found by experiments or theories. The 

equation (la) involves only nondimensional parameters and is therefore 
Independent of the plate parameters. Using equation (1), the nonlinear static 
and dynamic behaviors of a plate can be predicted and they are applicable to 
plates of any size, boundary conditions, material properties. 

From the plot of P versus q for static condition (fig. 1) regions of 
hardening and softening spring behavior are found and there are also regions 
of negative stiffness, (e.g. , between A and B). Notice that for R - 1, 

dynamic motion starting from C will pass through A and B and ends up 

at C’. 

Static equilibrium positions are found by putting p - 0 in equation (1) 
and correspond to the points where the curve crosses the q-axis in figure 1. 
Note that for R > 0, there are three equilibrium values of q; R 1 / 2 , 

— r 1/2 an( j 

zero. The last value, zero, is an unstable position as the ‘ 
stiffness is negative. We can rewrite equation (1) as 
3 2 

q - q Q q ■ p (2a) 



(2b) 


q/ft 2 + 2? q/ft + (q 3 - q*q) - p 
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where q Q ■ equilibrium position - R^/^. 

From equation (2a), dp/dq * 0 when q « *1/^3 q Q , these are the end 

points of region of negative stiffness, e.g., A,B in figure 1 


Undamped Free Vibration of Postbuckled Plate 

From equation (2b), neglecting damping effect and external force, the 
free vibration is obtained as 


1 »' 3 2 

-y q + (q - q Q q) - o 

ST o 


(3) 


With the substitution *q = q , equation 3 can be integrated. For the 

**r 

initial conditions q(o) * 0 ana q(o) * q (q is the initial position, q > q ) 

s s so 

at t * 0 the solution can be written 


•2 A f 2 2w 2 , 2 „ 2^ 

q * — lq s - q J lq s + q - 2 q Q J 


(A) 


2 2 2 

With the expression q + q - 2 q set equal to 1, the above formula is a 

s o 

standard linear vibration equation for initial value problem. 

The other extreme point of vibration, after the system was released at 


q s , q e , is found by substituting q = 0 in equation (A) to obtain 


2 2 1/2 

q - — q or q = (2 q* - qV /Z 

e s e os 


(5) 


When q < ^2q , 2q - q >0, 2q - q <q (since q > q ), for the 
n s o n o s ’ o s o s n o 

2 2 2 2 

second solution, 0 < q < q and q - q * q - q , thus the oscillation 
* e o o e s o 

Is around q , from q (> q ) to q (< q J . The first solution q * -q 
n o s o e o e s 

is on the negative side and thus not reached practically when q < 

s o 


A 



/— 2 2 

When q > >2 a , 2q - q <0, the second solution gives an imaginary 

s o o s 

number, thus the only possible solution is q e « -q s , which means that the 
oscillation is from q s to -q g . Also the motion passes through both 
equilibrium positions of q c and -q 0 and q - 0 is the new mean 
position, instead of -the original static value, q 0 . 

t— 2 2 

When q « *2 q , 2q -q *0, thus q - 0. 

From (4), by substituting q * it can be written as. 


/2 dq 



By numerical integration of (6) for the motion between the extreme points, 

q 8 , <le» time histories for the free vibration of various amplitudes were 

determined for a plate with q Q « 1 and are shown in figure 2. When 

q < , it can be seen that the period of vibration increases as amplitude 

s 

increases and the motion is not symmetrical about the static value, 1. When 
q - ^2 ■ 1.414, the period is theoretically infinity as it 

takes infinite time to approach zero. However, when q > /2 the period will 

s 

decrease with increase of amplitude. The displacement also passes through 

both equilibrium positions, q D * 1 and q 0 * -1 , thus indicating 

snap-through motion. The change of mean position is q c (from q D 

to 0). Also the rms value is found to be approximately 1 when q s *= 1.5. 

Essentially, the postbuckled plate shows softening spring behavior initially 

and after snap-through motion accompanied by a change of equilibrium position 

it shows a hardening spring behavior. 

The free vibration response characteristics reported in reference 2 also 

shows that resonance frequency decreases to zero when snap-through motion is 

initiated and the frequency increases when the magnitude of snap-through 
motion gets larger. , 



Random Vibration 


The method of equivalent linearization can be used to solve the nonlinear 
forced vibration equation with damping (from eq. (2b)); 

1 " C • t 3 2 \ 

q + 2 ^ q + (q - q Q q) - P (7) 

2 

The mean square displacement of a buckled plate, <q > due to white noise 
excitation with spectral density S pp can be obtained for small and large 
magnitudes, as described in Appendix II. 

4 , 

For small excitations a < 0.45 q [a = nfls /44, a is a nondimensional 

o PP 

force parameter) , there is no snap-through motion and 

<q 2 > - [~q 2 + (q* + 3ct) 1/2 ]/3 .and q - q Q . (8) 

4 

For a 2 q Q there is persistent snap-through motion in almost every 
cycle of oscillation and 

<q 2 > = [q 2 + (q* + 12 a) 1/2 ]/6, and q = 0 . (9) 

4 4 

For 0.45 q Q a < 2 q Q , snap-through motion is intermittent, the mean 
position as well as the mean square values are very unsteady. However, the 
mean square value can be taken approximately from interpolation between the 
two end points — the point of no snap-through and persistent snap-through 
motion. 

From equations (8) and (9) the variation of <q 2 > with the excitation 
parameter <* for different values of compression parameter X is shown in 
figure 3. The rate of increase of response with excitation is highest when 
intermittent snap-through motion starts. When persistent snap-through motion 
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is attained, the response increases much more slowly with increases in 
excitation, showing hardening spring behavior. 

From figure 3 and equations (8) and (9), the variation of response with 
compression parameter for various levels of excitation are plotted as shown in 
figure 4. For a given excitation level, the response increases with increases 
in compression as it approaches the buckling point. After initial post- 
buckling, persistent snap-through occurs and the response continues to 
increase until a certain point for which only intermittent snap-through motion 
can be induced. After that point, the response decreases with further 
increases in compression. The point of maximum displacement corresponds to 
the point for which the excitation is just sufficient for persistent 
snap-through motion. Also, the point of maximum response occurs at 
increasingly greater plate curvatures, or larger values of u/u c , as the 
excitation level increases. These trends have also been predicted in a 
qualitative description by Jacobson. 5 

EXPERIMENTAL RESULTS 
The Test Set-Up 

Dynamic tests on several 0.032-inch thick aluminum plates were carried 
out using point excitation at the center of the plate by an electromagnetic 
shaker. Displacement response was measured by strain gauges. As shown in 
figure 5, excitation force was applied in both directions but without applying 
any bending constraint on the plate using rounded point screws connected 
together to the shaker by a rectangular frame around the specimen. This 
direct attachment method of excitation was used instead of the base excitation 
of the supporting frame because it can ensure large excitation force in the 
low-frequency range (0-20 Hz). 
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Results for Discrete Frequency Excitation 

In a preliminary test, it was found that snap-through motion was most 
readily excited by an excitation frequency of 5 Hz. The variation of the 
rms value of strain response and reduction of the mean value of strain due to 
snap-through with excitation level are plotted in figure 6. There is a region 
of unsteady response as snap-through motion is initiated. When more power is 
put into the shaker, the excitation force decreases with increase of 
response. This is accompanied by a comparable reduction in the mean value of 
the oscillation approximately equal to the original static value, as predicted 
in figure 2. 

Random Responses 

Broadband (0—100 Hz) excitation was used to excite random response 
(fig. 7). The low-frequency response (0-20 Hz) dominates the strain response 
when snap-through motion is initiated at u/u c - 4.0 (fig. 7) and the 
fundamental modal response is not evident. The dominance of low-frequency 
response for large amplitudes of post-buckled plates was also reported in 
references 5 and 6. 

Figure 8 shows the transfer functions (strain/force) of a flat plate with 
no buckling load, u/u c - 0, for various excitation levels. The resonance 
peak broadens and the overall level decreases as excitation level is 
increased. This is a typical hardening spring behavior. 

Figure 9 shows transfer functions of a buckled plate. The resonance peak 
also broadens and decreases as did the flat plate, but the response at low 
frequency increases as excitation level is Increased. ' Increases in the low- 
frequency response are due to the onset of snap-through motion for which the 
natural frequency is very low (fig. 2). 
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The variation of mean square strain parameter (square of the ratio of 
dynamic strain to static strain at R * 1) with compression parameter is shown 
in figure 10. The general trend agrees well with predicted results from the - 
single-mode formula (fig. 4) and the points of maximum responses are near the 
curve for static values. However, there is a large discrepancy between the 
experimental results and theoretical prediction for an excitation level of 
<* ■ 6, which indicates that the single-mode representation used in the 
analysis over-predicts the stiffness of the buckled plate. More modes may be 
required to represent the deformation pattern and give a lower overall 
stiffness value. 

CONCLUSIONS 

A simple formula was derived for the transverse load-displacement 
relationship of a plate under in-plane compression and compared with results 
from experiment using electromagnetic excitation method on an aluminum plate 
with postbuckling deflection. The comparison shows that the simple formula 
predicts the general trend of the highly nonlinear behavior of snap-through 
motion under dynamic excitation. The general characteristics of snap-through 
motion are: 

(1) The mean position of the oscillation is zero; 

(2) The r.m.s displacement value is approximately equal to the static 
equilibrium value of the plate when snap-through motion is just 
initiated; 

(3) Snap-through motion is most readily excited by low-frequency 
excitation; and 

(4) For a given random excitation level, the maximum response is found in 
the postbuckling configuration when the excitation is just sufficient 
for persistent snap-through motion. 
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The identification of the nonlinear characteristics found in post-buckled 
plates should be very helpful in studying the corresponding characteristics in 
other curved plates such as cylindrical panels. 
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APPENDIX I 


Derivation of Single-Mode Formula 
For Any Plate Under Compression 

The starting point is the Von Karman equations of large deflection of 
plates. The tranverse equilibrium equation and the in-plane comptability 
equation can be expressed, respectively, as 


and 


where 



(A— 1 ) 


(A-2) 


F ■ Airy stress function 

°x» a y> T xy ” in-plane stresses 

w ■ transverse displacement 

D ■ plate flexural stiffness * Et 3 /12(l-v 2 ) 

E - modulus of elasticity 
t “ plate thickness 
v ■ Poisson’s ratio 

h ■ transverse load per unit area, i.e. the pressure. 

The steps are the same as in the Rayleigh-Ritz method described in reference 4 
and summarized as follows: 
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1 . Transverse Displacement: 


w(x,y) - tQ<Kx)4'(y) 

where Q is the modal displacement coefficient 

<Kx) ’Ky) is the buckling mode shape function. 

The corresponding modal force coefficient, P, is 
b a 

t / / h(x,y) <Kx) ’Ky) dxdy . 

0 0 

Stress Function : substituting w from (A-3) into (A-2), 

r ■ " 2 « 2 [ | f « V’> V y) - H ' I P y 


(A-3) 


where 


(A-4) 


4>^(x), ’l'j(y) are higher order functions related to <Kx) , 'l'(y), 


respectively, fjj depends on i,j (details in ref. 4); and 
P x , Py are mean compressive stresses in the x and y direction. 
Mean Compressive Stresses: for the edge displacement 




(A-5) 


From (A-4) and (A-5), 


P 

x 


where 


P 


u 


P 

u 




(A-6) 


P u is the mean compressive stress due to edge shortening if the plate 
is flat and C X y is a constant related to *Kx), ♦(y)« 

Similarly, with edge displacement in y direction being zero, 

P - vp - E C Q 2 (A-7 ) 

y x y 

Cy is a constant related to t(y). 
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4. In-plane Strain Energy : 


/ / -|=r la 2 + a 2 - 2v a a + 2(1 + v) t J dxdy 

j j 2E 1 x y xy xy J 

o o 


From (A-4), (A-6) and (A-8), 

V - fp 2 + P 2 - 2vP P ) + Etab eQ* 
e 2E ' x y x y J 

where e Is a constant related to fjj In (A-4). 

Bending Strain Energy ; 

Dt r b t a (3 2 w 3 2 w-v 2 . , 

v b‘— I > '■72 7* y 

o o 3x 3y 


(A-8) 


(A-9) 


(A-10) 


From (A-3) and (A-10), 


V. » DQ 2 d (A-ll) 

D 

where d is a constant related to $(x), tp(y). 

6. Linear Static Equilibrium Equation : 



This equation applies for flat plate with small value of Q, for which 
V e can be neglected. 

Substituting (A-9) into the above equation, 

2DdQ - P 

defining K ■ 2Dd, the linear modal equilibrium equation is obtained. 

KQ - P (A-12) 

K is thus the linear modal stiffness. 
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7. Nonlinear Static Equilibrium Equation : 


3(Ve + Vb) _ 

3Q 

Substituting from steps 4 and 5 and rearranging. 



Etab n 2 
Dd 4 


P tabC 
u 

Dd 


^ + 1 ] 2Dd = P 


Dd 

Substituting ^ P c 

xy 

2Dd - K, K is the linear modal stiffness from (A-12). 



(A-13) 


( A-14) 


the equation (A-14) becomes 

Q - (/->)] K-P (A-15) 

Q P 

2 2 

This is the nonlinear modal equilibrium equation. Q /Qp * s tbe effect 

of large displacement, P u /P c is tbe effect of compression. 

u P 

Defining 2 a “ P c ’ and 'F' “ u“ ’ substitutin 8 <1 = Q/Qp» P “ p / K Q p > 

(l - v ) a c c 

The general nondimensional equation is obtained by dividing (A— 15) by 
KQ p . 

q^ -(u/u c - 1 ) q * p (A-16) 

Putting p=0 for loading with in-plane compression only, 

P 

For jr- < 1 < 1), q “ Q “ 0 (A-l 7) 
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For 


T y K > »• ^2 ■ <r’ '> 

C c Qp C 


or 


(— - 1 ) 


(A-18) 


From the above it can be seen that P<» and u c are the critical mean 
compressive stresses and edge shortening respectively, thus when 

Q « or q * 1 (A-19) 


- 2 


8 . 


Therefore, Qp is the value of Q when u/u c - 2 
Kinetic Energy 


1 b ® ,o 

V T * J pt / / w dx dy 

o o 


(A-20) 


where p is the density. 
Substituting (A-3) to (A-20) 


v-H 2 


O b a O 2 

where M * pt / / <J> (x) i|> (y) dx dy M is the modal mass, 
o o 


(A-21) 


9. Lagrangian Equation for Linear Vibration 


^b + i_ 

3Q 3t ^ 


(A-22) 


Using results from (A-12), (A-20), the equation of motion is obtained, 

KQ + MQ* * P CA— 23) 

The natural linear circular frequency of the flat configuration is thus 
given by 


ST = K/M 


( A-24) 
15 



The modal mass M can be written in terms of ft 


a 


(A-25) 


10. Lagrangian Equation for Nonlinear Vibration 

3(V + V. ) . 3V 

— ® L_ + — — — = p 

3Q 3t 9Q 


(A-26) 


using results from (A-15) and (A-23) the modal equation of motion is 
obtained , 


[ ~ “ • ~ 1 ] KQ + MQ* = P 

QP c 


(A-27) 


Dividing both sides by KQp and using results from (A-16) and (A-25), the 
nondimens ional equation of motion is obtained. 


3 0 . 1 

q - Rq + —j q 

a 


p (A-28) 

The damping effect can be similarly included by dividing the modal damping 


force 2McaQ ( £ is the modal damping coefficient) by KQp, which gives 
2 1 • 

— q. Thus the nondimensional equation of motion can be written as 


+ la + ,3 


q - Rq 


(A-29) 


a 
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APPENDIX II 


Derivation of Nonlinear Acoustic Response Using 
Equivalent Linearization 

Flat Plate 

The linear equation of motion for the flat plate is 


1 x, • 

“2 q + 2-jjq + q- p 
2 

The mean square value of q, <q >, due to white noise excitation with 

irAS 

spectral density Spp is given by <q Z > ■ — w ^ ere “ is a 


(B-l) 


nondimens ional excitation parameter. 

Postbuckled Plate without Snap-Through Motion 

The nonlinear equation of motion for the postbuckled plate is 

1 •• t • 3 2 

-o q' + 2 ^ q + (q - q q) - p (B-2) 

Q 2 11 ° 

Substituting q * q 0 + Aq into equation (B-2), (Aq is the dynamic 
displacement around the static value q Q ) , the equation become 

K A^*+ 2 £ A<5 + (Aq 3 + 3Aq 2 q + 2A qq 2 ) = p (B-3) 

QiC o o 

The equivalent linear equation is 

Aq + 2 Aq + kAq « p (B-4) 

Si L u 

3(Aq^ + 3Aq 2 q + 2Aqq 2 ) 

where k *= — — 



Thus from eq. (B-4) 


<V> 


irflS 


_E£ 


4Ck 


3<4q 2 > + 2q 2 Q 


(B-5 ) 


Solving eq. (B-5) for <Aq > 


<Aq > - [-q + C q n + 3°) J /3 


(B-6) 


Post buckled Plate with Snap-Through Motion 

The mean position of q, is zero, thus equation (B— 2) is used directly and 
the equivalent linear equation is 


1 •• C • 

-=-q + 2-jjq + kq-p 


(B-7) 


where k 

From eq. (B-7) 


3 (q 3 - q 3 q) 

5S 


2 . 2 
3<q > - q o 


irftS 

« 2 > ■ -dF 


3 <q 2 > - <£ 


(B-8) 


Solving eq. (B-8) for <q > 

<q 2 > - [q 2 + (q* + 12 ^ /2 h 

The Excitation Level for Snap-Through Motion 


(B-9) 


2 2 

For (B-9) to be valid, <q > q Q so that most of the oscillation consist 

of complete snap-through motion (the r.m.s. of which is approximately 

2 2 2 
q Q , as shown in fig. 2). By substituting <q > from (B-9) to <q > >^ q 0 


a > 2 q Q is obtained as the condition for persistent snap-through motion. 
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If <q 2 > < 1/3 q 2 , the oscillation rarely gets into the region of 
negative stiffness (between A and B in fig. 1), there is no snap-through 
motion. By substituting <q 2 > from (B-6) to <q 2 > < 1/3 q 2 , we obtain 

a < 0.45 q is obtained as the condition for no snap-through motion. 
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Figure 1. Plots for transverse force equilibrium equation of 
a plate under in-plane compression. 
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Figure 3. The variation of displacement response with excitation level in a postbuckled 
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COMPRESSION PARAMETER (u/u 




plate specimen 



Figure 5. The schematic diagram for the dynamic test of a plate 






strain 
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Figure 6. The variation of strain response with excitation level in a plate. 
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Figure 7. The strain response to random excitation in a postbuckled plate with 
persistent snap-through. 
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Figure 8. The transfer function of the strain response of an isotropic plate 
under random excitation (no compression load u/u = 0.0). 



aluminum alloy 4"x8"x.032" BUCKLED PLATE 



eojoj. jo Ajisuep |ejjoeds 
/eo-ioj pue ujB4S jo Ajjsuep jBjpsds ssojq 


random excitation (u/u 



aluminum alloy 4"x8 w x.032 



J0J0UJBJBd U|B4S 0JBnbS UB0^ 


29 


Figure 10. Variation of random strain response with increasing compression in an 
isotropic plate. 
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